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Abstract 

Possible dark states could be induced after derivations of the entrainment of matter 
induced by a surface wave propagating along the flexible vacuum-matter boundary by con- 
sidering the nonlinear coupling between the interface and the rarefaction effect. The non- 
relativistic limit of the relativistic Navier-Stokes equations was considered and analytically 
solved by a perturbation approach. The critical reflux values associated with the product 
of the second-order body forcing and the Reynolds number (representing the viscous dissi- 
pations) decrease as the Knudsen number (representing the rarefaction measure) increases 
from zero to 0.1. We obtained the critical bounds for possible dark states corresponding to 
specific Reynolds numbers (ratio of wave inertia and viscous dissipation effects) and wave 
numbers which might be linked to the dissipative evolution of certain large-scale structure 
during the relativistic heavy-ion collisions. 
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1 Introduction 

Recently a new state of matter has been created in Au-|-Au collisions at RHIC, and surpris- 
ingly it was found to flow as a perfect fluid [1]. The kinematic shear viscosity of this nearly 
perfect fluid has been determined and found to be rather small compared to conventional low- 
temperature fluids. Relevant researches in heavy-ion physics has focused on constraining the 
transport properties of hot and dense nuclear matter using experimental data from RHIC and 
studying relativistic hydrodynamics for viscous fluids in order to describe the expansion of the 
fireballs created in relativistic heavy- ion collisions. 

Meanwhile the mean cosmic density of dark matter (plus baryons) is now pinned down to be 
only ca. 30% of the so-called critical density corresponding to a 'flat'-Universe. However, other 
recent evidence — microwave background anisotropics, complemented by data on distant super- 
novae — reveals that our Universe actually is 'flat', but that its dominant ingredient (ca. 70% of 
the total mass energy) is something quite unexpected: 'dark energy' pervading all space, with 
negative pressure. We do know that this material is very dark and that it dominates the internal 
kinematics, clustering properties and motions of galactic systems. Dark matter is commonly as- 
sociated to weakly interacting particles (WIMPs), and can be described as a fluid with vanishing 
pressure. It plays a crucial role in the formation and evolution of structure in the universe and 
it is unlikely that galaxies could have formed without its presence [2]. Analysis of cosmological 
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mixed dark matter models in spatially flat Friedmann Universe with zero A term have been 
presented before. A large majority of dark energy models describes dark energy in terms of the 
equation of state (EOS) pd = uj pd (cf. Refs. 3 and 4), where to is the parameter of the EOS, 
while pd and pd denote the pressure and the energy density of dark energy, respectively. The 
value C(j = — 1 is characteristic of the cosmological constant, while the dynamical models of dark 
energy generally have u; > — 1. The case of the growing cosmological term A and its implications 
for the asymptotic expansion of the universe and the destiny of the bound systems have been 
studied in Ref. 4 using above system of equations. Their results showed that even for very slow 
growth of A (which satisfies all the conditions on the variation of Gjv), in the distant future the 
gravitationally bound systems become unbound, while the non-gravitationally bound systems 
remain bound. 

Influential only over the largest of scales-the cosmological horizon-is the outermost species of 
invisible matter: the vacuum energy (also known by such names as dark energy, quintessence, 
x-matter, the zero-point field, and the cosmological constant A) (cf. Refs. 5 and 6). If there 
is no exchange of energy between vacuum and matter components, the requirement of general 
covariance implies the time dependence of the gravitational constant G. Thus, it is interest- 
ing to look at the interacting behavior between the vacuum (energy) and the matter from the 
macroscopic point of view. One related issue, say, is about the dissipative matter of the flat 
Universe immersed in vacua [7] and the other one is the macroscopic Casimir effect with the 
deformed boundaries [8]. 

Theoretical (using the Boltzmann equation) and experimental studies of interphase nonlocal 
transport phenomena which appear as a result of a different type of nonequilibrium representing 
propagation of a surface clastic wave have been performed since late 1980s (cf. Refs. 9 and 
10). These are relevant to rarefied gases (RG) flowing along deformablc elastic slabs with the 
dominated parameter being the Knudsen number (Kn = mean-free-path/L^, mean-free-path 
(mfp) is the mean free path of the gas, Ld is proportional to the distance between two slabs) 
[11-13]. The role of the Knudsen number is similar to that of the Navier slip parameter Ng (cf. 
Ref. 14); here, Ng = jiS/d is the dimensionless Navier slip parameter; S is a proportionality 
constant as Ug = St, t : the shear stress of the bulk velocity; Ug : the dimensional slip velocity; 
for a no-slip case, S = 0, but for a no-stress condition. S = oc, fi is the fluid viscosity, d is one 
half of the distance between upper and lower slabs). 

Note that, there are some models, like the MIT bag model and its descendants, where mat- 
ter is in a bag, in which there is no vacuum and the vacuum is outside. In this particular 
model one might speak of a clear vacuum-matter boundary. Here, borrowing the idea of the 
MIT bag model, the transport driven by the wavy elastic vacuum-matter boundary will be 
presented. The flat-Universe is presumed and the corresponding matter is immersed in vacua 
with the interface being flat-plane like. Wc adopt the macroscopic or hydrodynamical approach 
and simplify the original system of equations (related to the momentum and mass transport) 
to one single higher-order quasi-linear partial differential equation in terms of the unknown 
stream function. We then introduce the perturbation technique so that we can solve the related 
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boundary value problem approximately. To consider the originally quiescent gas for simplicity, 
due to the difficulty in solving a fourth-order quasi-linear complex ordinary differential equation 
(when the wavy boundary condition are imposed), we can finally get an analytically perturbed 
solution and calculate those physical quantities we have interests, like, time-averaged transport 
or entrainment, perturbed velocity functions, critical unit body forcing corresponding to the 
possible dark states. These results might be closely linked to the vacuum-matter interactions 
(say, macroscopic Casimir effects) and the evolution of the Universe (as mentioned above : the 
critical density [2]). Our results also show that for certain time-averaged evolution of the matter 
(the maximum speed of the matter (gas) appears at the center-line) there might be existence of 
negative-pressure states. 

2 Formulations 

The matter is presumed to be a fluid associated with a shear viscosity but no bulk viscosity 
and no heat conduction here (the geometrized units arc adopted so that G = c = 1 and the 
Einstein's field equations : G/^i, = STrr^y). The stress tensor (for this fluid) is [15] 

Here, poi ^o, P, and are the rest-mass density, specific internal energy, pressure, and the fiuid 
4-velocity, respectively, r] is the coefficient of viscosity and is related to the kinematic viscosity 
u hy r] = pqu. a^i, is the shear tensor (the detailed expression could be traced in Ref. 15). In 
general a F-law equation of state P = {T — l)poeo could be presumed. Thus, we have basic fluid 
variables 

p* = poau^e^'^, e* = (poeo)^^^ au^e^'^ , Sk = p*huk, 

where is the the conformal exponent and h = 1 + eo + P/po is the speciflc enthalpy (a is 
the elapse or proper time elapsed in moving between the neighbouring spatial hypersurfaces). 
The conservation of stress-energy [15] gives TJ^'^ = and the law of baryon number conservation 
V^(/9ow'*) = gives the relativistic continuity, energy, and Navier-Stokes equations 

dtp* + di{p,v') = 0, (1) 
de, + diie^v') = ^ae'''^riipoeo)^^-^^/^a'"^aap, (2) 

dtSk + diiSkv') = -ae^-^Pfc + 2{ae'"^na'^),^ + ae'^'^gfivaap - ^pohuaUp), (3) 

where f * = u^/uq is the 3- velocity. The quantity is determined by the normalization condition 
u'^Ui, = 1, which has 

w' = pI + e-"^rSiSj[l + ]-\ 

P:t^{we'^'P / p^y 
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with w = p^av9 {'jij is the spatial or 3- metric). We remind the readers that the stress tensor 
T'"^ generates the source terms in the field evolution equations : 

Sij = -^^i^3 + P^ij - ^Wij, 

where /?* is the shift (for gauge conditions) or displacement in spatial coordinates in moving 
between the neighbouring spatial hyper surf aces. 

Equations (1-3) arc too difficult to be solved analytically or even by a perturbation approach. In 
this work we only consider the nonrelativistic limit of above equations. Meanwhile the vacuum 
is presumed to be incompressible (cf., e.g., Ref. 16, the speed of sound propagating in this 
'vacuum' is formally very large (or infinite) rather than zero as in the empty vacuum; this 
implies the incompressible vacuum). 

Note that the first theories of relativistic dissipative fluid dynamics are due to Eckart [17] and 
to Landau and Lifshitz [18]. The difference in formal appearance stems from different choices 
for the definition of the hydrodynamical four-velocity. These conventional theories of dissipative 
fiuid dynamics are based on the assumption that the entropy four-current contains terms up 
to linear order in dissipative quantities and hence they are referred to as first order theories 
of dissipative fiuids. The resulting equations for the dissipative fluxes are linearly related to 
the thermodynamic forces, and the resulting equations of motion are parabolic in structure, 
from which we get the Fourier-Navicr- Stokes equations. They have the undesirable feature that 
causality may not be satisfied. That is, they may propagate viscous and thermal signals with 
speeds exceeding that of light. Extended theories of dissipative fluids due to Grad [19], Miiller 
[20], and Israel and Stewart [21] were introduced to remedy some of these undesirable features. 
These causal theories arc based on the assumption that the entropy four-current should include 
terms quadratic in the dissipative fluxes and hence they are referred to as second order theories 
of dissipative fluids. The resulting equations for the dissipative fluxes are hyperbolic and they 
lead to causal propagation of signals [22] . A qualitative study of relativistic dissipative fluids for 
applications to relativistic heavy ions collisions has been done using these first order theories. 
The application of second order theories to nuclear collisions has just begun [23]. 
The flat-plane boundaries of this matter-region or the vacuum-matter boundaries are rather 
flexible and presumed to be elastic, on which are imposed traveling sinusoidal waves of small 
amplitude a (possibly due to vacuum fluctuations). The vertical displacements of the upper 
and lower interfaces {y = h and —h) are thus presumed to be r] and — ry, respectively, where 
r] = acos[27r(x — ct)/X], A is the wave length, and c the wave speed, x and y are Cartesian 
coordinates, with x measured in the direction of wave propagation and y measured in the 
direction normal to the mean position of the vacuum-matter interfaces. The schematic plot of 
above features is shown in Fig. 1. 
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It would be convenient to simplify these equations by introducing dimensionless variables. We 
have a characteristic velocity c and three characteristic lengths a, A, and h. The following 
variables based on c and h could thus be introduced : x' = x/h, y' = y/h, u' = u/c, v' = v/c, 
rl = r]/h, ip' = ^/(c/i), t' = ct/h, p' = p/{pc?), where ip is the dimensional stream function, u 
and V are the velocities along the .x- and y-directions; p is the density, p (its gradient) is related 
to the (unit) body forcing. The primes could be dropped in the following. The amplitude ratio 
e, the wave number a, and the Reynolds number (ratio of wave inertia and viscous dissipation 
effects) Re are defined by 

a 2'Kh ch 

e = -, a = — -, Re = — . 
h \ V 

We shall seek a solution in the form of a series in the parameter e : = ipo + eV'i + + ■ ■ 

dp/dx = {dp/dx)o + e{dp/dx)i + e^{dp/dx)2 + with u = dijj/dy, v = —dip/dx. The 

two-dimensional (x- and y-) momentum equations and the equation of continuity could be 

in terms of the stream function ip if the p-term (the specific body force density, assumed to 

be conservative and hence expressed as the gradient of a time-independent potential energy 

function) is eliminated. The final governing equation is 

and subscripts indicate the partial differentiation. Thus, we have 
d 1 

— V^V'o + V'Oi/V^Vox - 'ipOx^'^i^Oy = '^^'^^0, (5) 
d 1 

— v^V'i + V'oj/V^Vix + V'lj/V^Vox - V'oiV^V'ii/ - V'liV^Voi/ = (6) 

V'Oa;VV22/ - -^Ix^^^^ly - '^2x'^^'^0y = ;^^^^2, (7) 

and other higher order terms. The (matter) gas is subjected to boundary conditions imposed by 
the symmetric motion of the vacuum-matter interfaces and the non-zero slip velocity : n = =F 
Kn du/dy (cf. Refs. 11,12, and 13), v = ±dr]/dt at y = ±(1 + r/), here Kn=mfp/(2/i). The 
boundary conditions may be expanded in powers of r] and then e : 

V'o?/|i + e[cosa(x-t)'0O2/?/|i +V'i3/|i] cos^ q;(x - t) -|- V'2y|i+ 

cosa(x - t)tpiyy\i] -\ = -Knlipoyyli + e[cosa(x - t)tpoyyy\i +ipiyy\i\+ 

^2^ i'0yyyyh ^^^2 _ ^ ^^^^^^ _ + ^^j/yli] + •••}, (8) 

V'Oxli + e[cosQ(x - t)V'oxy|i + V'lxli] + e'^[ ^^'''^^^^ cos^a(x - t)+ 

cosa{x — t)ipixy\i + '(p2x\i\ + • ■ ■ = —ea sin a(a; — t). (9) 
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Equations above, together with the condition of symmetry and a uniform {dp/dx)o, yield : 

V^o = i^o[(l + 2Kn)y-^], Ko = ^{-^)o, (10) 

^1 = + <^*(y)e-^"(^-*)}, (11) 

where the asterisk denotes the complex conjugate. A substitution of ipi into Eq. (6) yields 

(P (f 
{— r - + iaRell - KqU - + 2Kn)]}(— r - a^U - 2iaKoRe(j) = 

or if originally the (matter) gas is quiescent : Ko = (this corresponds to a free (vacuum) 
pumping case) 

(P dp 

i-r^ - o?)i-r^ - oPU = 6? = oP-iaRe. (12) 

The boundary conditions are 

0j^(±l) ±^j/y(±l)Kn = 2i^o(l ±Kn) = 0, 0(±1) = ±1. (13) 

Similarly, with 



we have 



Dyyyy = -^-^^^^ly - ^*^vy)y^ (15) 

[— . - (4^2 - 2iaRe)]i—^ - Aa^)E - i2aReKo(l - + 2Kn) 
dy^ dy^ 

(-^ - Aa^)E + iAaKoReE + '■^^{(t)y<l)yy - (p(pyyy) = 0; (16) 



and the boundary conditions 

Dyi±l) + ^[0,,(±1) + ryy{±l)] " 2i^o = TKn{ 1 [<^,,, (±1) + .^^,,(±1)] + Dyy{±l)}, (17) 

^2/(±l) + \<l>yyi±i) - ^ = TKn[^</),,,(±l) + ^,,(±1)], (18) 

E{±l) + ^<Py{±l) = (19) 

where Kq is zero in Eqns. (16-18). After lengthy algebraic manipulations, we obtain = coe"^-|- 
cie-"2/ + C2e^f + C3e-"f , where co = {A + Ao)/Det, a = -{B + BQ)/Det, C2 = (C + CQ)/Det, 
C3 = -(T + To)/Det; Det = Ae" - Be"" + Ce" - Te"", 

^ = e°^a2(r2e-2" - s^e^^) - 2aae-''w + aae°^z(e-2"r + e^^s), 



7 



Bo = e°^a2(r2e-2« - s'^e'^^) - 2aae-'^w + aae'^ z{e-'^^r + e^^s), 
C = eT'^aoiiwse^-'^ - rze'^'^) - ae'^"+'^z{az - as) + ae-'^w{ae^-'^w - ae'^'^r), 
Co = e'^aaiwse'^-'^ - r^e""") - az{zae^''-^ - ae"s) + aw{ae~^'^+'^'^^w - ae-(2"+")r), 
T = e-'^aa{zse^+'' - rwe-^''+^^) - aa{e'^'^-^zr - e^ws) + a'^e-''{-e'^'^z'^ + e'^'^w'^), 
To = e'^aa{zse'^+" - ru;e-("+")) - aaie'^zr - e'^'^'^ws) + a'^e'^{-e'^'^z'^ + e'^^'w'^), 
with r = (1 — aKn), s = (1 + aKn), w = {1 — aKn), z = (1 + aKn). 

To obtain a simple solution which relates to the mean transport so long as only terms of O(e^) 
are concerned, we see that if every term in the x-momentum equation is averaged over an interval 
of time equal to the period of oscillation, we obtain for our solution as given by above equations 
the time-averaged (unit) body forcing 

where XIq is the integration constant for the integration of equation (15) and could be fixed 
indirectly in the coming equation below. Now, from Eq. (17), we have 

Dy{±l) ± KnDyyi±l) = " ^ [</>,,( i 1 ) + .^^.(il)] T Kn{i[</>,,,(±1) + </>;,,(±l)]}, (21) 
where Dy{y) = Hoy^ + aiy + 02 + C{y), and together from equation (15), we obtain 
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with gi = a + a* , g2 = a + a, gs = a — a* , g4^ = a — iy, = a — a* , gQ = a + a* . In 
realistic applications we must determine ITo from considerations of conditions at the ends of the 
matter-region, oi equals to zero because of the symmetry of boundary conditions. 
Once Ho is specified, our solution for the mean speed {u averaged over time) of matter-flow is 

C/ = = ^{C(y) - C(l) + i?o - KnCj,(l) + Uo[y^ - (1 + 2Kn)]} (23) 

where Rq = —{[(pyyi^) + 0^j/(l)] ~^^[4'yyy{^) +</'yj/j/(l)]}/2, which has a numerical value about 
3 for a wide range of a and Re (playing the role of viscous dissipations) when Kn= 0. To 
illustrate our results clearly, we adopt U (Y) = u{y) for the time-averaged results with y = 1" in 
the following. 
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3 Results and Discussion 

We check our approach firstly by examining Rq with that of no-shp (Kn= 0) approach. This 
can be done easily once we consider terms of Dy{y) and C{y) because to evaluate Rq we shall at 
most take into account the higher derivatives of ^(y), like <f>yy{y), 4>yyy{y) instead of 4>y{y) and 
escape from the prescribing of 02- 

Our numerical calculations confirm that the mean streamwise velocity distribution (averaged 
over time) due to the induced motion by the wavy elastic vacuum-matter interface in the case of 
free (vacuum) pumping is dominated by Rq (or Kn) and the parabolic distribution — no(l — y^). 
Rq which defines the boundary value of Dy has its origin in the y-gradient of the first-order 
streamwise velocity distribution, as can be seen in Eq. (17). 

In addition to the terms mentioned above, there is a perturbation term which varies across the 
channel : C{y) — C{\). Let us define it to be 

^^^) = ;^[^(^^-^^')] ^''^ 

We remind the readers that the Reynolds number here is based on the wave speed. The physical 
trend herein is also the same as those reported in Refs. 12 and 13 for the slip-fiow effects. The 
slip produces decoupling with the inertia of the wavy interface. 

Now, let us define a critical refiux condition as one for which the mean velocity U{Y) equals to 
zero at the center-line y = (cf. Fig. 2). With equations (15,23-24), we have 



dp [a'Re'F{Q)/2m + KnC'{l)-Ro] 
= ^"fe^^ = -(l + 2Kn) ^^^^ 

which means the critical refiux condition is reached when IIq has above value. Pumping against a 
positive (unit) body forcing greater than the critical value would result in a backward transport 
(reflux) in the central region of the stream. This critical value depends on a, i?e, and Kn. There 
will be no reflux if the (unit) body forcing or pressure gradient is smaller than this IIq. Thus, 
for some Ho values less than IIo^^, the matter (flow) will keep moving or evolving forward. On 
the contrary, parts of the matter (flow) will move or evolve backward if Ho > IIo^^ . This result 
could be similar to that in Ref. 17 using different approach or qualitatively related to that of 
Ref. 4 : even for very slow growth of A, the gravitationally bound systems become unbound 
while the nongravitationally bound systems remain bound for certain parameters defined in Ref. 
4 (e.g., 77). 

We present some of the values of no(Q;, i?e;Kn = 0,0.1) corresponding to dark states which 
satisfy j\ U{Y)dY = in Table 1 where the wave number (a) has the range between 0.20 and 
0.80; the Reynolds number {Re)= 0.1, 1, 10, 50, 100. Here, the dark state is directly related to 
the zero-mass-flux (in average) state in bulk sense although there are locally velocity differences 
across the (confined) space. It means the gravity effects are not absent locally considering the 
velocity distribution across the confined space. To be specific, the dark state is difficult to be 
observed since there is no mass flux even though this state happens! 

We observe that, from Table 1, as Kn increases from zero to 0.1, the critical Hq or time-averaged 
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(unit) body forcing decreases significantly. For the same Kn, once Re is larger than 10, critical 
reflux values Ho drop rapidly and the wave-modulation effect (due to a) appears. The latter 
observation might be interpreted as the strong coupling between the vacuum-matter boundary 
and the inertia of the streaming matter-flow. The illustration of the velocity fields for those dark 
states arc shown in Figure 2. There are three wave numbers : a = 0.2,0.5,0.8. The Reynolds 
number (Re) is 50. Both no-slip and slip (Kn=0.1) cases are presented in Fig. 2. The line of 
value U = in Fig. 2 is schematic and could represent the direction of positive and negative 
velocity fields. 

Some remarks could be made about these dark states (or solitons) : the matter or universe being 
freezed in the time-averaged sense for specific dissipations (in terms of Reynolds number which 
is the ratio of wave-inertia and viscous effects) and wave numbers (due to the wavy vacuum- 
interface or vacuum fluctuations) for either no-slip and slip cases. This particular result might 
also be related to a changing cosmological term (growing or decaying slowly) or the critical 
density mentioned in Ref. 2. If we treat the (unit) body forcing as the pressure gradient, 
then for the same transport direction (say, positive x-direction) , the negative pressure (either 
downsdtream or upstream) will, at least, occur once the timc-averaged flow (the maximum speed 
of the matter (gas) appears at the center-line) is moving forward! 

Meanwhile, the time-averaged transport induced by the wavy interface is proportional to the 
square of the amplitude ratio (although the small amplitude waves being presumed), as can be 
seen in Eqn. (12) or (20), which is qualitatively the same as that presented in Ref. 9 for analogous 
interfacial problems. In brief summary, the entrained transport (pattern, cither postive or 
negative and there being possible dark states) due to the wavy vacuum-matter boundary is 
mainly tuned by the (unit) body forcing or Ho for fixed Re (viscous dissipation). Meanwhile, 
nOcr depends strongly on the Knudsen number (Kn, a rarefaction measure) instead of Re or a 
(wave number). We hope that in the future we can investigate other issues [22-26] using the 
present or more advanced approach. 
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Table 1: Dark states values (Hq) for a flat vacuum-matter boundary. 







Re 










Kn 


a 


0.1 


1 


10 


50 


100 





0.2 


4.5269 


4.5269 


4.5231 


4.4496 


4.3275 


0.5 


4.6586 


4.6584 


4.6359 


4.4086 


4.2682 


0.8 


4.9238 


4.9234 


4.8708 


4.5714 


4.4488 


0.1 


0.2 


2.4003 


2.4000 


2.3774 


1.9532 


1.2217 


0.5 


2.4149 


2.4132 


2.2731 


0.7728 


-0.9054 


0.8 


2.4422 


2.4379 


2.0718 


-0.5885 


-3.4151 
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Fig. 1 Schematic diagram of the deformable motion of the vacua-matter boundary. 
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Kn=0.0, a=0.8 
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Kn=0.1. a=0.8 
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Fig. 2 Demonstration of the dark states : the mean velocity field U (Y) for 

wave numbers a = 0.2, 0.5, 0.8. The Reynolds number is 50. Kn is the rarefaction measure 

(the mean free path of the particles divided by the characteristic length) . 

The U = line is schematic and illustrates the directions of positive and negative U{Y). 

The integration of U(Y) w.r.t. Y for these velocity fields gives zero volume (mass) flow rate. 



